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•  A  longer  service  life  of  the  structure  operating  in  extreme 
environments  consisting  of  severe  thermal  and  mechanical  loading 
conditions. 

•  Combine  the  properties  of  two  dissimilar  Materials 

•  The  load  carrying  capacity  of  plate-type  structure  is  enhanced. 

•  An  increased  operational  life  of  the  structure. 
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2.  Basic  Assumptions  and 
Preliminaries 


Figure  1.  A  simply  Supported  functionally  graded  (FG) 
plate  shown  in  2-D,  under  an  explosive  blast 
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2.  Basic  Assumptions  and 
Preliminaries-Continued 
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Basic  Assumptions 

•  The  plate  mid  surface  is  referred  to  a  Cartesian  Orthogonal  System  of 
Coordinates  (x,  y  z). 

•  Zis  the  thickness  coordinate  measured  positive  in  the  upwards  direction 
from  the  mid-surface  of  the  plate  with  h  being  the  uniform  plate  thickness. 

•  The  plate  is  assumed  to  be  thin.  As  a  result,  the  Kirchoff  Assumptions, 
of  zero  transverse  shear  stresses,  will  be  assumed  leading  to  the 
classical  plate  theory. 

•  The  functionally  graded  plate  is  composed  of  ceramic  and  metal  phases 
whose  properties  vary  smoothly  and  continuously  across  the  wall  thickness 
from  one  surface  to  another. 
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Based  on  a  rule  of  mixtures,  the  following  form  of  the  variation  of 
mechanical  properties  (Young’s  Modulus,  Density,  and  Poisson’s  Ratio) 
are  postulated  as 


Denotes  a  Gene 


^(z)  -  (^ceramic "  ^metal)K;(z) +  ^metal 


ic  Property 


Volume  Fraction 


Specific  Values  of  the  Respective  Properties 


Note :  For  Vc(z)  =  0,  P(z)  =  Pm  and  For  Vc{z)  =  1,  P(z)  =  Pc 

=^c(z)^[0,1] 
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Type  1 :  Symmetric  Distribution  of  the  Constituent  Phases 


The  phases  vary  symmetrically  through  the  wall  thickness  in  the  sense 
of  having  full  ceramic  at  the  outer  surfaces  of  the  plate  and  tending 
toward  full  metal  at  the  mid-surface 


Where, 


sgn(z)  = 


V(z,  N) = 

0,  ifz  =  0 
-1 ,  if  z  <  0 
+  1,  if  z>  0 


/  \N  .  ,  \N 

z  |  1+sgnz  z  1-sgnz 


yh/2, 


+ 


V 


-h/2 


J 


Signum  Function 


N  =  Volume  Fraction  Index  (0  <  N  <  oo) 


UNCLASSIFIED:  Dist  A.  Approved  for  public  release 

5  August  2010 


GVSETS 


tHTJim  oi 


M 

3.  Types  of  FGM  Models  -Continued 

?  *  MODELING  AND  SIMUimiON,  TESTING  and  VRLIOHTION 


V 


9 


Pc  (Ceramic  Phase)  Pm  (Metallic  Phase) 

~  V  A  h/2 

A  ▼ 

Figure  2.  Distribution  of  the  constituent  materials  through 
The  plate  thickness  for  the  symmetric  case. 

at  z  =  +  h/2,  for  any  N,  V  =  1,  and  as  a  result  P(  ±  h/2 )  =>  Pc,  and  for  z = 0, 

V  =  Oand  P(0)  =3  Pm. 

Type  2:  Asymmetric  Distribution  of  the  Constituent  Phases 

The  phases  vary  non-symmetrically  through  the  wall  thickness  and  in 
this  case  there  is  Full  ceramic  at  the  outer  surface  of  the  plate  wall  and 
full  metal  at  its  inner  surface. 
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3.  Types  of  FGM  Models-Continued 
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V(z,  N)  =  [(/7  -  2z)/2h]N 


(A-l) 

Pc  (Ceramic  Phase)  Px  (Metallic  Phase) 


i 


hjl 


-h/2 


Figure  3.  Distribution  of  the  constituent  materials  through 
The  plate  thickness  for  the  asymmetric  case. 


at  z=  +  h/2,  P(h/2)  =>  Pm,  and  for  z=  h/2,  P(h/2)=>  Pc.  At  the  Midsurface, 
z  =  0,  and  for  k  =  \  P(0)  =  ^(PC  +  Pm) 
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3-D  Displacement  Field 


u(x,y,z,t)  =  u0{x,y,t)-z 

v(x,y,z,t)  =  v0(x,y,t)-z 
w(x,y,z,t)  =  w0(x,y,  t) 


8w0(x,y,t ) 
dx 

dw0(x,y,  t) 
dy 


Uq,  Vq,  Wq  =>  2-D  displacement  quantities  of  the  mid-surface 
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Non-linear  Strain-Displacement  Relationships 


SXX  - 


du  1  (  dw \ 

- 1 —  - 

dx  2  v  dx  J 

dv 


Syy  - 


1  (  dw 

H —  - 

dy  2\dy 


J 


Yxy 


£xz 


i  ii  dv  dw  dw 

- 1 - 1 - 

dy  dx  dx  dy 

£yz  —  £ZZ  ~  0 
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4.  Theoretical  Developments- 
Continued  ha  iv 
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Substitution  of  the  displacement  field  into  the  strain-displacement  relationships 
Gives  the  3-D  strain  measures  across  the  plate  thickness  as 
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4.  Theoretical  Developments- 
Continued 
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Plane  Stress  Constitutive  Equations 


G  XX 

On 

0)2 

0 

£xx 

Gyy 

<,  — 

0[2 

Q22 

0  • 

£  yy  > 

Gxy, 

0 

0 

°66_ 

* 

j 

Gxz  ~  Gyz  ~  Gzz  ~  0 


Where, 


Oil  -  ^22  -  “  2  ’  ^66  -  ZZZT.  L’^16’^26  ~  0 


1-v‘ 


1  -v‘ 


2(1  +  v) 
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Governing  Equations 


fj  ( SU  +  SV-  SK)dt  =  0  =>  Hamilton' s  Principle 

*0 


t0,  ^  are  two  arbitrary  instants  in  time 

U,  denotes  the  strain  energy 

V,  denotes  the  work  done  by  surface  tractions,  edge  loads,  and 

body  forces 

K,  denotes  the  kinetic  energy  of  the  3  -  D  body  of  the  Structure 
S,  denotes  the  variational  operator 
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Strain  Energy 


•hi  2 


5U  =  -b0  Uh/  2^xxSsxx  +  ayySeyy  +  axysrXy  )dzdnO 


External  Work 


sv  =  - 


Jq0  Pt(X’y)^X,y,^dn0-lx^/2{^*yySv  +  ^*yxSu  )dzdx 


\ 


h/2 

y!-h/2 


ivi 


* 


* 


V 


c 7xxSu  +  cr XySv  | dzdy 


Kinetic  Energy 

■h/2 


SK  =  \QJ^2p(z)WSWdzdCl0 


5  August  2010 


16 


GVSETS 


VlfTJtm  ai 


■ 


sv  4.  Theoretical  Developments- 
Continued 
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Equations  of  motion  in  terms  of  Stress  resultants 


dNxx  oN 


dx 


+ 


xy  


dy 


=  0 


dNww  dN 


•yy 


dy 

d2M 


+ 


xy  


dx 


=  0 


XX 


+  2 


82M 


xy 


d2M 


+ 


yy 


dx 

d  f 


dxdy  dy1 


+ 


d 

dx 


f 


V 


m  di/i/n  .. 

N yy - 1"  N y\/ 

xx  dx  xy 


dw0 ^ 

8y  , 


+ 


dy 


Nww  — ^  +  N 


dl/l/0^ 


V 


'yy 


dy 


xy 


dx 


+  Pt  -  Cw0  =  IqWq 


J 
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4.  Theoretical  Developments- 
Continued 


MODELING  rnq  SIMULATION,  TESTING  and  VALIDATION 


Where, 


( Njj .  Mjj )  -  j  h^2 n-jj (1.  z) c/z  (i,j  =  x,y,xy)  =>  Stress  Resultants  and 

Stress  Couples 

Simply  Supported  Boundary  Conditions 


wQ  =  Mxx  =  Nxy  =  0,  Nxx  =  Nxx  on  x  =  0,  L, 

WQ  =  Myy  =  NyX  =  0,  N yy  =  Nyy  Ony  =  0,  4 

For  Compressive  Edge  Loading,  N*xx  =  -N%x  and  N*yy  =  -Nyy 
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The  first  two  equations  of  motion  are  satisfied  by  assuming  the  following 


d2<j> 

dxdy 


It  can  be  shown,  the  details  of  which  are  not  shown  here,  that  the  third 
equation  of  motion  can  be  expressed  in  terms  of  displacements  and  a 
stress  potential  as 


DV4i/i/0 


r d2(j)d2\N§  ^  d2l/l/0 

vdx2  dy2  dxdy  dxdy 


d2(j)  d2W0^ 
dy2  dx2  y 


+  l0w0  +  Cw0  = 
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5.  Solution  Methodology- 
Continued 


Where, 


C-£1£3-Ej 

£,(1-v2) 

One  more  equation  in  terms  of  the  transversal  deflection  and  the 
stress  potential  is  needed.  This  will  come  from  the  strain  displace¬ 
ment  relationships,  by  eliminating  the  in-plane  displacements,  which 
is  known  as  the  compatibility  equation  which  can  be  shown  to  be 
given  by 


v V  =  5 


-2  A 

d w0 

dxdy 


d2wQ  d2wQ 
dx2  8y2 
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-j2  ^2 

Where,  V2  =  — ^  + 


dx2  dy2 


Laplacian  Operator 


Assume  a  solution  form  for  </>,  Wq 


w0(x,y,t)  =  wmn(t)s\nAmxs\nMny 


4  —  Amn(t)  cos  2/Lmx  +  Bmn(t)cos2/jny  +  Cmn(f)cos2Amxcos2//ny 
+  Dmn(t)  sin  2/tmxsin  2  juny  + \  N*xxy2  + \  Nyyx2 
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Where, 


/trp  =  iun  =  nnjL-z,  m,n  =  1, 2, 3 etc. 


Utilizing  the  compatibility  equation  and  comparing  coefficients  gives 


2  2  2  2 

a  ( +\  -  ^wmn(^)Mn  d  n\  _  n  /+\  _  n  u\  -  n 

™mnv)  ~  o  ’  & mnV )  ~  _ _  9  » ~  ^mnV)  ~  ^ 


32/l‘ 


m 


32ju 


n 
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Continued 
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Transversal  Pressure 

Let, 

P?  =  ^mn(0  sin /tmxsin//ny 
By  integrating  over  the  plate  area  gives 


Pmn(t) 

or 


4  /  / 

7T Jo2#  Pf (0 sin Amxsin /jnydxdy 
l1l2 


16Pf(t)  .... 

— Where,  m,n  =1 

n 
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'  5.  Solution  Methodology- 
Continued 


MODELING  and  SIMULATION,  TESTING  and  VALIDATION 


Finally,  applying  the  Galerkin  method  to  the  third  equation  of  motion 
Gives  the  nonlinear  equation  differential  equation  governing  the 
Structural  response  of  FG  plates  under  external  excitation. 


Wmn(t)  +  2Amn<ymnl/i/mn(f)  +  <z>mnwmn(t)  +  ^ mnwmn( 0  —  Pmn(t) 


wmn(t)  =$  Amplitude  of  Deflection 

a)mn  =  -I^rnnUo  =>  The  natural  frequency  of  the  plate 

Amn  =  C/2l0a)mn  =>  Nondimensi  onal  damping  factor 

^ mn  ~  E-\  (2m  +  fin  )/l  6 
Pmn(f)  =  16Pf(f)//0^-2 
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'  5.  Solution  Methodology- 
Continued 
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4  A  2 

UTT  /  4  0  2  2  2  4  4\  /Vw^  ,  2  2^  2\ 

Kmn  =  — j-{m  +2 mrn  y/  +  n  y/  )  +  —^, — (m  +  n <by/  ) 

L\  ^-1 

Where, 

y/  =  mt-2  =>  Aspect  Ratio 

O  =  Nyy  / A/xx  =>  Compressive/Tensile  Edge  Load  Ratio 


The  governing  nonlinear  differential  equation  is  solved  using  the 
4th-Order  Runge-Kutta  method  assuming  zero  initial  conditions. 
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6.  Blast  Loading 


For  a  free  in-air  spherical  air  burst,  the  pressure  profile  over  time  is  given 
in  figure  4  as 


Fig  4.  Incident  Profile  of  a  blast  wave 
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6.  Blast  Loading-Continued 
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\\ 


The  wave  form  shown  in  figure  4  is  given  by  an  expression  known  as 
The  Friedlander  equation  and  is  give  as 


Pt(t)  =  (Pso  ~  P0) 


f  .  .  \  -a 

“I  _  f  ~  *3 


V 


t 


t-t  ^ 

1  ‘a 

lp  j 


P  J 


Where, 


Pso  =  =>  Peak  Overpressure  over  ambient 


’3  72  Z 


Ft  is  the  Standoff  Distance 

Z=R/ WV3  scaled  distanced  W  is  the  equivalent  weight  of  charge 

of  TNT  in  terms  of  kilograms 
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6.  Blast  Loading-Continued 

MODELING  m  SIMULATION,  TESTING  and  VALIDATION 

UNCLASSIFIED:  Dist  A.  Approved  for  public  release 


P0  is  the  ambient  pressure 
ta  is  the  time  of  arrival 

tp  is  the  positive  phase  duration  of  the  blast  wave 
t  is  the  time 


For  conditions  of  STP  at  sea  level,  the  time  of  arrival  and  the  positive 
phase  duration  can  be  determined  from 


Arrival  time  or  positive  phase  duration 

1 


W}3 


'■Mj 


Cube  root  scaling 


Arrival  time  or  positive  phase  duration  for  a  reference  explosion  of  charge  weight,  l/l/-, 

It  should  be  noted  that  the  standoff  distances  are  themselves  scaled 
According  to  the  cube  root  law 

GVSETS 
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Validation  of  the  Present  Approach 

Comparisons  are  made  (Akay,  H.U.,  1980)  who  considered  explosive  blast  of 
Metallic  plates  under  a  step  loading  excitation  given  by 

Pt(t)  =  P[H(t)-H(t-t0)] 

H(t)  is  referred  to  as  the  Heaviside  Step  Function  defined  as  H(t)  =  1  for  t  >  0 
andH(t)  =  Ofor  t  <  0 

The  geometrical  and  material  properties  used  were 

L\  =  2.438/77,  h  =  0.00635/77,  y/  =  1 

Em  =  70.3 GPa,  pm  =  2547kg/  m3,  vave  =  0.25 
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Fig  5.  Comparisons  of  solutions  of  the  time  history  of 
the  Central  deflection  Under  a  Step  Load  P  =  48.82Pa, 
t$  =  0.2  sec 
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Present  Results 

A  ceramic-metal  functionally  graded  plate  with  the  following  properties  were 
used 

Table  1 .  Material  Properties 


Metal 

(TS-6AI-4V) 

Ceramic 

(Aluminum 

Oxide) 

E,  Modulus  (GPa) 

105.7 

320.24 

u,  Poisson’s  Ratio  (Unitless) 

0.2981 

0.26 

p,  Density  (Kg/m3) 

4429 

3750 

m  =  1m,  iff  =  1 ,  h  =  0.0254  m) 
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Table  2.  Air  blast  parameters  versus  distance  for  a  one 
Kilogram  TNT  spherical  air  burst  (Kingery  and  Bulmash, 
1984). 


Standoff 

Arrival 

Positive  Phase 

Distance, 

Time, 

Duration, 

(m) 

(msec) 

(msec) 

1.0 

0.532 

1.79 
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7.  Results-Continued 
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N= 0  (Fully  Ceramic) 
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Fig  6.  Implications  of  the  volume  fraction  index  on  the 
deflection-time  history  of  an  asymmetric  FG  Plate  without 
damping. 
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Fig  7.  Implications  of  the  symmetry  on  the  deflection-time 
response  of  FG  plate.  =0.1 
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Fig  8.  Implications  of  various  amounts  of  damping  on  the 
deflection-time  response  of  an  asymmetric  FG  Plate.  (N=0.5) 
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Fig.  9.  The  implications  of  the  compressive/tensile  edge  loading 
on  the  deflection-time  response  of  an  asymmetric  FG  plate. 

(Au  =0.2,  A/ =  0.5) 
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Fig.  10.  The  effects  of  various  aspect  ratios  of  an  asymmetric  FG 
plate  on  the  deflection-time  response.  (A/=0.5) 
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Fig.  11 .  The  implications  of  the  value  of  the  decay  parameter 
on  the  time-deflection  response  of  a  asymmetric  FG  Plate 
without  damping.  (N=0.5) 
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•  A  Structural  Model  of  the  dynamic  response  of  functionally  graded  plates  exposed 
to  a  free  in-air  blast  with  simply  supported  boundary  conditions  has  been  presented. 

•  As  the  volume  fraction  of  ceramic  increases,  the  deflections  become  smaller  with 
higher  frequencies  in  contrast  to  metal  phase. 

•  The  symmetric  construction  due  to  being  stiffer  has  smaller  deflections  and  higher 
frequencies  as  compared  to  the  asymmetric  construction 

•  For  a  fixed  compressive  edge  load,  as  the  metal  concentration  increases,  the 
deflections  increase  with  decreasing  frequencies. 
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